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Abstract 

We study Lagrangian trajectories and scalar transport statistics in decaying Burgers turbulence. We 
choose velocity fields solutions of the inviscid Burgers equation whose probability distributions are specified 
by Kida's statistics. They are time-correlated, not time-reversal invariant and not Gaussian. We discuss 
in some details the effect of shocks on trajectories and transport equations. We derive the inviscid limit of 
these equations using a formalism of operators localized on shocks. We compute the probability distribution 
functions of the trajectories although they do not define Markov processes. As physically expected, these 
trajectories are statistically well-defined but collapse with probability one at infinite time. We point out 
that the advected scalars enjoy inverse energy cascades. We also make a few comments on the connection 
between our computations and persistence problems. 
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1 Introduction 



Lagrangian trajectories driven by a velocity field u(x, t) are solutions of the differential equation: 

dx(t) 



dt 



u(x(t),t) (f) 



As known from the precursors, Richardson, Kolmogorov, Batchelor,... [Q, they acquire peculiar properties 
when the flow becomes turbulent. These properties are probably going to play an important role in the 
understanding of fully developed turbulence. For example, the recent proof Q of the existence and uniqueness 
of the stationary state for the inviscid forced Burgers turbulence is based on an analysis of these trajectories. 

Statistical properties of these trajectories may be deciphered by looking at transport phenomena in tur- 
bulent systems. Recent studies of the Kraichnan's advection models ||] have made these expected properties 
more explicit. Kraichnan models assume that the velocity fields is Gaussian and white-noise in time. These 
simplifications lead to the solvability of the models. See refs.|| for recent studies of the Kraichnan models 
for incompressible fluids and refs.^j, |E| for compressible ones. Two kinds of behavior have been observed: 

(1) Statistical ill-definedness, meaning that two trajectories starting at the same point have a non vanishing 
probability to be far apart at later time, 

(2) Trajectory collapse for compressible enough fluids, meaning that two trajectories starting initially at 
different positions have a non-zero probability to follow the same path after some time, 

However: 

(3) Properties (1) and (2) do not seem to occur simultaneoulsy. 

The motivation of the present work is to decipher whether these properties arc more robust and hold true 
for more realistic velocity fields than those chosen in Kraichnan's models. Of course we could not solve the 
problem with a velocity field describing a real three dimensional turbulent system. Instead we shall consider 
(unrealistic) velocity fields, solutions of the Burgers equation which in 1+1 dimensions takes the form: 

d t u + ud x u — vd 2 x u = (2) 

where u = u(x, t) is the (compressible) velocity field and v the viscosity. This is a variant of the Navier-Stokes 
equation in which the role of the pressure has been neglected. Although we shall stick to one dimensional 
space, some of the following considerations could be generalized to higher dimensions. 

No external force is applied to equation (Q) . So its inviscid limit v — ► corresponds to decaying turbulence 
whose statistical description consists in finding probability distribution of the velocity fields solution of eq. (|^) 
given random initial data. One usually expects a more universal behavior at large time. Thus, we shall 
consider a family of velocity fields, solutions of the inviscid limit v — > of the Burgers equation, whose 
probability distribution describes the long time behavior of large classes of initial conditions. These velocity 
statistics are those first introduced by Kida ||. In contrast to the Kraichnan model, the velocity fields are 
then not white-noise in time, not time-reversal invariant and not Gaussian. 

For compressible fluids, one may look at two kinds of transport phenomena depending on whether one is 
looking at the advection of a tracer, that we shall denote by T(x, t), or at the advection of the density of a 
pollutant, that we shall denote by p(x, t). The corresponding viscosity is written n. The equations governing 
these transports are: 

d t T(x, t) + u(x, t)d x T{x, t) - nd 2 x T(x, t) = (3) 
d t p(x,t) +d x (u(x, t)p(x,t)) ~ ndlp(x,t) = (4) 

They differ by the order of the derivative and velocity. 

In the inviscid limit v — ► 0, solutions of the Burgers equation develop shocks at which the velocity is not 
smooth. This non-smoothness implies that the naive definition of the trajectories does not apply. Therefore, 
these trajectories and the transport equations have to be dealt with carefully. As we shall see, a correct 
definition of the transport equations will turn out to be: 

d t T(x,t) + ^(u(x+,t)+u(x-,t)^d x T(x,t) - nd 2 x T(x,t) = (5) 

d t p(x,t) + d x ^(u{x+,t) + u(x-,t) S )p(x,t) S ) - Kdlp[x,t) = (6) 



with u(x ± ,t) — lim e ^o+ u i x i e i0- Although equations (^,g) seem to be naively equivalent to equations 
), they are not since in the inviscid limit the velocity field u(x,t) is not smooth. 
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In the limit k — > 0, eqs.(g,g) have a natural interpretation in terms of Lagrangian trajectories. However, 
the naive equation (|l|), which is actually meaningless since u(x, t) is discontinuous, has to be modified into: 

^l\ + = ±(u(x(t)+,t)+u(x(t)-,t)) (7) 

Again this differs from eq.(Q) because u(x,t) is not smooth. The physical meaning of this modification is 
clear. At points of discontinuity, the ill-defined velocity is replaced by the velocity of the shock which, as is 
well known since long ago [0, ||, is just the average of the velocities just before and just after the shock. 
Once this will be done, we shall describe how to compute the probability distribution functions (p.d.f.) of 
the trajectories and we shall use them to discuss the properties of the transport equation (|5|) in the limit 

K — ► 0. 

This paper is organized as follows. In the following section we recall basic facts concerning the Burgers 
equation and the velocity profiles we shall use. Section 3 is devoted to give a precise the definition of 
Lagrangian trajectories in the inviscid limit v — > and to the relation with the correct form of transport 
equations and their solutions. In Section 4 we establish identities, called equations of motion, which are 
valid inside correlation functions. This is based on operators localized at shocks and their algebra. In 
Section 5, the backward and forward probability distribution functions of the trajectories are introduced 
and their formal properties emphasized. We use the identities established in Section 4 to verify that these 
p.d.f. are solutions of the transport equations. In Section 6, we make explicit computations for one and two 
particle distributions. We check the consistency with the expected physical properties of the trajectories. 
In particular, we show by different approaches that the trajectories are statistically well-defined but that 
particles have a non-vanishing probability to collapse. This is in agreement with the general properties of 
Lagrangian trajectories mentioned above as (1),(2) and (3): As we deal with a highly compressible fluid, the 
alternative (2) is realized and the alternative (1) is excluded. The connexion with persistence problems is 
made. Finally arguments indicating that the energy cascade in scalar advection in these flows is inverse, i.e. 
towards the large scale, are presented in Section 7. 



2 Velocity profiles 



This short section is devoted to the specification of the statistics of the velocity profiles to be used in this 
paper. 

• In order to fix notations, we recall a few elementary facts concerning Burgers equation (see e.g. |t], |) 
and II and references therein). As is well known, the equation is solved by implementing the Cole-Hopf 
transformation which maps it to the heat equation. This works as follows. Let Z(x,t) = exp[— ^j$(x, t)} 
where t) — d x ${x, t). Eq. (||) for u is mapped into the heat equation for Z: 

[d t - vdl] Z{x, t) = 0. 

Thus, given the initial condition u(x, t = 0) = uq(x), the velocity field at a later time t is recovered from the 
potential &{x,t) given by the relation 



exp 



2v 



dy 



exp 



(8) 



with <&o(x) standing for the initial potential such that uo(x) = d x &o(x). The inviscid Burgers equation 
corresponds to the limit v — » 0. The solution is then given by solving a minimalization problem: 



u(x, t) = d x $(x, t) with t) = min $ (?/) + 

y 



- yf 
2t 



(9) 



Outside shocks the minimum is reached for only one value ?/» ofy, the solution of the equation uo(y*)t — x—y*. 
The velocity is u(x,t) — x ~?* — uo(y*)- It is effectively a local solution of the inviscid Burgers equation 
since, by the minimum condition defining y*, we have u(x,t) — uq(x — tu(x,t)). A simple geometrical 
construction of the solution (^) is described in refs.[^, ||. For large t, coincides approximately with one of 
the local minima of &o(y) and it practically does not change under small variations of x so that, in between 
the shocks, the velocity is approximately linear with the slope j. 
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Shocks appear when the minimum is reached for two values y\ and y 2 of y. Let $12 = $0(2/1,2) be the 
value of the initial potential at these points. Then eq. (|^) allows one to determine the velocity profile u s (x, t) 
around and inside the shocks at finite value of the viscosity v by expressing exp [— j^$ s (x, t)] as the sum of 
contributions from the two minima. One obtains: 



u s {x,t) = ~ [x- +2/2M - 7^ tanh (-^ (x-£ a t - + y 2 ) 



(10) 



where ^ s = yi - y 2 > and £ s - '''' ' l ' 
^ — * 0, cq.(n3) becomes 



. The width of the shock is of order L ~ — . In the inviscid limit 

V1-V2 c Ms 



u s (x,t)l =0 = £ S T^+ X t for ±(x-x s (t)) >0 



(11) 



where a; s (t) = !; s t + 5(2/1 + J/2) is the time £ position of the shock which moves with the velocity £ s and 
follows a Lagrangian trajectory. The values of the velocity on the two sides of the shock are: 



Ms 

2t 



(12) 



so that tj- is the amplitude of the shock. 

• To mimic this large-time behavior, following Kida M, we choose as velocity profiles the ansatz: 

[x-V]) 2 ^ 



u(x, t) = d x S(x, t) with S(x, t) — min ( <pj + 



2t 



(13) 



The points {<t>j,yj)jez specify a given realization. For any realization, i.e. for any data of the points 
(cj)j,yj), these ansatze ( |l3| ) are solutions of the inviscid Burgers equation. They have exact sawtooth profiles^ 
with slope 1/t. In this ansatz all shocks are created at time t — 0. The later time evolution is then governed 
by the shock collisions. Thus different times are strongly correlated. 

Following Kida we shall concentrate on velocity statistics specified by demanding that (<f>j, yj)jez be 
a Poisson point processQwith intensity J = e^dipdy. 

This choice of statistics ensures that the velocity u(x,t) is self-similar with characteristic length l(t) ~ \ft 
which means that su{sx,s 2 t) = u(x,t). Here and in the following, = means an equality in law, i.e. inside 
any correlation functions. We could as well choose other intensities for the Poisson process. This amounts 
to choose other scalings for the characteristic length. 




The sawtooth velocity profile. 



1 In particular, the velocity is not defined by the above formulae at the shocks. This is at the origin of most of the following 
discussions. 

2 The basic rules to manipulate such processes are briefly recalled in appendix ^ and [b| where some explicit computations 
are made. 
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3 Lagrangian trajectories and transport equations. 



Lagrangian trajectories x(t) starting at point xq at time to are denned as solutions of the evolution equation: 

dx(t) 



dt 



— u(x(t),t) with x(to) — xq (14) 



In this section, we shall specify the equation governing Lagrangian trajectories in the inviscid case. This 
requires first a detour through the v ^ situation. 

• The above differential equation is well-posed for a velocity field u(x, t) solution of the Burgers equation 
with finite non vanishing viscosity v ^ 0, since then u{x,t) is smooth enough. 

However the limit v — ► is delicate: 

— If the point x(t) of a trajectory is far from shocks, the velocity is then regular around that point even 
in the inviscid limit and the trajectory is well defined. At large time, the velocity far from shocks is of the 
form u{x,t) = j(x — y„) with approximately constant and the trajectories are then straight lines. This 
applies as long as the trajectories are away from shocks. 

— Assuming that shocks are diluted, the trajectories near a shock in the inviscid limit v — > may be 
analyzed using the velocity profile (fy]). In this environment, solutions of the Lagrange equation x = u(x, i) 
are such that: 

sinh (—^-(x(t) — x s (t))j exp ( — — s - ] = constant. 
V Aft J V ovt J 

with x s {t) the time t position of the center of the shock. Recall that the width of the shock is of order 
l c ~ — . This equation means that particles away from the shock take a finite time to enter the shock. Once 
they are in the shock they move coherently with it with velocity almost equal to £ s = x s (t) . But they never 
cross the shock center. 

• If we want to recover the v = limit behavior directly in the inviscid case with the ansatz (|l3|) for 
the velocities, we have to be careful. At discontinuities of the velocity, eq.(|l4|) does not make sense for two 
reasons : the velocity is not defined at the shocks and the derivative of a differentiable function cannot 
exhibit pure discontinuities. 

A simple modification that will ensure the gluing of particles to shocks, the main feature at finite but 
small viscosity, is the following: 

— First we define u(x,t) = ^(u(x + ,t) + u(x~,t)). For the ansatz (|f3|), this definition makes sense for any 
x and extends the definition of u(x,t) to shocks (obviously u = u away from shocks). 

— Then we demand that trajectories be continuous and satisfy 

dx(t) ^ X (t + e)-x(t) = _ { ^ t) 



dt ' e^0+ 



c 



If we assume that the shocks form a discrete set (no limit points)^ these two requirements ensure that 
trajectories are uniquely defined for t > to once the boundary condition x(io) = £o is specified. Since the 
velocity of a shock is the mean of the velocities at the points just preceding and just following it, eq.([l5|) 
ensures that particles stick to shocks. 

• According to the ansatz (|T^), away from shocks, u(x, i) — j(x — y) for some y. So the trajectory is: 

x(t) = xq + (t — to)— — -, away from shocks (16) 

with xo the position at time to. This is true up to the time at which the particle meets a shock. Shocks 
are at the points where two parabolae 0i,2 + ^ X ~2t' 2 ^ minimizing eq.fjfj]) intersect. They move with a 
velocity £12 = Zyl ■ I n the time interval during which the shock exists, the trajectory equation is the shock 

x ( t ) = \{V~L + 2/2) + 62 t, on the shock (17) 

Once a particle is on a shock it follows it and the cascade of shocks arising from its collisions. Note that 
when two shocks hit they merge into a third shock. In particular, a particle not on a shock at time t has 

3 The probability distribution for the velocities ensures that this happens with probability one. 



_ 01—02 

12 

equation: 
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never met a shock before^ A general feature of the trajectories is that particles move at constant velocity 
on intervals of the form [t, t'[ with t < t' . 

• This definition of the Lagrangian trajectories ensures the physical fact that the velocity field in the 
inviscid limit is transported by the fluid. Indeed, since a particle moving along a Lagrangian trajectory keeps 
its velocity for a finite time interval e with e sufficiently small, one has: 



lim - 



u(x + eu(x, t),t + e) — u(x, t) 







Accordingly, the transport equation for a tracer T(x, t) moving in the inviscid velocity field u{x 1 i) will be: 

i r 



lim 

e-»0+ £ 



T(x + eu(x, t), t + e)- T(x, t) 



= 



(18) 



It coincides with the k — > limit of eq. (J5|) provided the Lagrangian trajectories are specified as in (|15[). 

We show now how the above equation ([l8]) can be solved. The idea is to find an implicit formula for the 
Lagrangian trajectories, taking any number of shocks into account. Fix xq and to, and consider the function 
X(x, t) = x — xq — (t — to)u(x,t) for t > to and x arbitrary. It is readily checked that ■*) > 

so that for fixed t, X(x,t) is a strictly increasing function of x with \mi x ^± 00 X(x, t) — ±oo. This means 
that we can define a function x(t) for t > to by the condition that X(x(t) + , t) > > X(x(t)~, t). It is 
cumbersome but straightforward to check that x(t) is the solution of ( pL 5|) with initial condition x(to) = Xo. 
Hence the solution of ( |l8|) with initial condition T(x,to) = 0{x — xq), where 9(x) is the Heaviside step 
function, is T(x, t) = 6{x — xo — (t — to)u(x, t)). By linearity, the solution with initial data T(x, to) = To(x) 
is T(x,t) — To(x — (t — to)u(x,t)). This solution develops discontinuities^ at the shocks, even if the initial 
condition is smooth. 



4 Operator localized on shocks and equations of motion. 

In this section, we discuss what happens to the Burgers equation in the inviscid limit. We shall argue that 
the actual inviscid Burgers equation is not the naive limit v — ► of eq.(0) but is: 



d t u(x, t) + - (u{x + ,t) +u(x ,t)^j {d x u(x, t)) 



= (19) 



with the equality valid inside correlation functions with velocity fields (with or without derivatives) away 
from x and velocity fields (without derivatives) at the point x. This is not quite the usual way to write the 
inviscid Burgers equation. So we shall start with the more familiar formulae and show the equivalence with 
(p9|). The argument will be based on an analysis of operators localized on shocks which may be used to 
derive equations of motion valid inside any correlation functions. 
• At v 7^ the Burgers equation (||) could be written as: 

(dt + ud x - v d 2 x + A 2 v {d x uf\ e A " = (20) 

Since c Au is finite in the inviscid regime, the distribution d x e Xu is well defined in this limit, and vd x e Xu 
vanishes when v — ► 0. Eq ( |20| ) can be rewritten in this limit as: 

(d t u[x,t) + u{x,t)d x u{x,t)^je Xu{x ^ + Ae(i,t)e Al, ^» S (21) 

Here e(x, t) is the dissipation field defined by e(x, t) = lim v(d x u) 2 . The product u(x, t)d x u(x, t) is ill-defined 
since it is a product of distributions. Eq.([2lj) should actually be read as: 

[d t + Xd X jd x ^e Xu ^^ + A 2 e{x, t)e x < x ^ = (22) 

This is the well-known inviscid equation of motion. The fact that the dissipation field survives the inviscid 
limit is sometimes called the dissipative anomaly. 

4 A different proof of the same result can be found in appendix |b] where it appears as a natural part of the argument. 
5 But no nastier singularities. 
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• The presence of shocks is at the origin of universal features which are independent of the details of the 
statistics. As explained in rcf.|l(|, they may be analyzed by looking at fields localized on the shocks. By 
definition, these fields may be represented for any realization as: 

O g (x,t)= 5(x-x s (t)) (23) 

shocks 

where the sum is over the shocks with x s (t) denoting the position of the shock, £ s its velocity and its 
amplitude. The function g(£, s , Ms) which specifies O g will be called the form factor of the operator. 

By using the velocity profile ( |l 0|) inside and around the shocks, one may map fields defined in terms of the 
velocity u(x, t) into the shock representation. The two basic examples described in ref. Jl0| are the generating 
functional (d x — j) e Xu ^ x ^ and e(x, t) e Xu i x ^) with e(x, t) the dissipation field. These fields are localized on 
the shocks. Indeed, outside shocks (d x u(x,t) — j) vanishes since away from shocks, u(x,t) — x ~ t v * with j/„ 
almost independent of x. Similarly, the dissipation field e(x), which is naively zero due to the prefactor v in 
its definition, is actually a non-trivial field since (d x u) 2 is singular in the inviscid limit. These singularities 
are localized on shocks and so is the dissipation field. In other words, dissipation takes place only at shocks. 

The shock velocity profiles (|l0|) at finite v can be used to regularize the ill-defined expressions that arise 
in a naive v = limit. In practice, given a local functionnal of the velocity, which is well defined at finite 
viscosity, one takes the v = limit in the distributional sense. For the above two examples, one obtains 
@: 

d x j) e^^) = -2]Te^ sinh(^i) S(x - x s (t)) . (24) 



and 



e(x,t)e^ = 2A- 3 £e^ cosh(^i) - smh(^)) S( 



x-x s (t)). (25) 



Now one may use the representation of the dissipation field as an operator localized on shocks to find 
alternative representations of them. Indeed eq.(p4[) implies: 

u(<u-i)e*«) 4) --|E eACs ( A ^ sinh ^) + ^r cosh ^)- sinh (^r)) *(*-*.(*)) 

However, looking at the product of the operator (|24|) with velocity at nearby points gives: 

±(u(x+,t) + u(x-,t)) f(d x u-- t )e x A =-|^e A ^A6sinh(^) S(x - x s (t)) (26) 

This is found using the fact that the velocity on the two sides of the shocks are uf — £ s T ■ Comparing 
these expressions with the form factor of the dissipation fields, eq. (p5|) , gives: 

Ae(a;,<)e AM(:E '* ) S Q (u(x+ , t) + u(x~ , <)) - u(x, tU (d x u(x, t)) e Xu( - x ^ (27) 

This is an extension of the well-known formula e(x) = lime?; [u(x) — u(x + I)] 3 . In this formula, ud x ue Xu 

means d\jd x e Au . As expected the dissipation field is located on the discontinuity of the velocity field. This 
relation is valid inside any correlation functions with other fields away from point x. 

Comparison of eq. ( p0[ ) with eq. ( p7| ) yields an alternative way of writing the inviscid Burgers equation in 
which the dissipation has completely disappeared: 



d t u(x, t) + ^ (u(x+,t) + u(x , i)) (d x u(x, t)) 



,A u{x,t) ^ ( 28 ) 



This is equivalent to eq.(^9[). It has a simple interpretation: it is the simplest possible point splitting 
regularization of the naive inviscid Burgers equation. The validity of this formula can be checked by hand 
in simple explicit correlation functions. 
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• Correlation functions of the velocity fields, without any derivative, are continuous as functions of 
the positions of the velocities. But the non-smoothness of the velocities in the inviscid limit implies that 
correlation functions of derivatives of the velocity field may be discontinuous and/or singular when points 
coincide. 

This has echoes on the products of operators localized on shocks: 

— Products of an operator localized on shocks times powers of the velocity field are discontinuous at coin- 
ciding points. These properties were illustrated in eq. (p6|) ; 

— Products of operators localized on shocks are singular at coinciding points. More precisely, fields localized 
on shocks form a closed algebra JlO| : 

O f (x,t) -O g (y,t) = S(x~y) O fg (x,t) + regular. (29) 

The contact term S(x — y) in this operator product expansion arises from the coinciding shocks in the double 
sum representing the product operator. This operator product expansion implicitly assumes that shocks are 
diluted. 



5 Lagrangian trajectory statistics. 

For non time-reversal invariant velocity fields one may consider backward and forward Lagrangian statistics. 

• The backward statistics encodes the probability distribution of the initial positions of the trajectories 
at time to knowing their positions at later time t > to. For n-trajectories they are given by the expectation 
values, 

n 

PlllXxj,t\x j ,to) = (Y[Vret.{Xj,t\x° j ,to)) with Vret. (x, t\x° , to) = 5{x° - x(t \x, t)) (30) 
3=1 

with x(to\x, t) the position of the trajectory at time to which will be at x at later time t > to- 

Although backward statistics are clearly a well-defined object from a probabilistic point of view, our 
representation of backward statistics, involving x(to\x, t) may seem inappropriate because trajectories may 
merge with increasing time, so that in general trajectories cannot be followed for decreasing time. However, 
the measure of the set for which one or more of the points Xj lies exactly on a shock at time t is zero, so 
that the backward trajectory is defined with probability one. This is obviously true as long as the points Xj 
are all distinct. When two or more of them coincide, things are not so clear. However, we shall be able to 
check explicitly that our backward statistics are well-normalized i.e. that: 

J ]Jdx j P^lXx j ,t\x j ,to) = l. 

3 

This ensures that we have not missed delta functions at coincident points. 

• The forward statistics codes the probability distribution of the final positions of the trajectories at 
time t knowing their initial positions at a previous time to < t. For n-trajectories they are given by the 
expectation values, 

n 

P lti X j' t \ X P t o) = (Y['Padv.{x j ,t\x j ,to)) With Vadv.{x,t\x°,t Q ) = 6{ X - x{t\x° ,t )) (31) 
3=1 

with x(t\xo,to) the position of the trajectory at time t which was at xq at the initial time to < t. Hence, 
x(to\x, i) and x(t\xo,to) are formally inverse functions: x(t\x{to\x,t),to) — x. The forward probability 
distribution functions are normalized such that: 

/ 11^ p il>3,*K°,*o) = i 

J 3 

• To deal with functions and not distributions, it is convenient to compute expectation values of products 
of step functions: 

n 

H^(x j ,t\x j ,to)=(Y[n{x j ,t\x j ,t )) with n(x,t\x°,t ) = 6(x° -x(to\x,t)) (32) 

3=1 
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with 9(z) the step function: 8{z) = for z < and 9(z) = 1 for z > 0. The functions give the 
probabilities for particles at points x 3 at time t to be at positions above x® at time to- They are such that: 



Pl n l{ Xj ,t\xlt Q ) = \{d x oH^\x 3l t\x%t ) 

3 



(33) 
(34) 



Since W n '(xj,t\Xj,ta) are expectation values of local functional of the velocity field not involving derivatives 
they can be computed directly from the velocity distribution functions. These are recalled in Appendix A. 

[nl 

Remark that P^. et will be regular at coinciding points since they do not involve derivatives of u, whereas 
^adv wn ^ ^ e sm S u l ar since they involve such derivatives. 

• Let us now argue that the backward statistics are related to the joint laws of the speeds u(xj,t), at 
least as long as the configuration is non-degenerate (no two points Xj coincide). In this case indeed, with 
probability one, no Xj lies on a shock, so each has a speed described by a single parabola, and then the 
same was true at any previous time. Hence with probability one, the particle passing at Xj at time t was at 
Xj — (t-to)u(xj, t) at time to (remember that as long as they do not meet a shock, particles move at constant 
speed). So only a trivial change of variables is needed to go from the joint law of the initial positions x® to 

Xi—X°- 

the joint law of the speeds u(xj,t), the relation being u(xj,t) — -j—f-f-- This ensures that the total mass 
of the backward distribution for non coincident points is unity, so that no finite probability is carried by 
degenerate configurations. This implies that 

H(x,t\x° ,t ) = 9(x° -x+ (t-t )u(x,t)). 

Hence the backward probability distribution is: 

V ret .(x,t\x°,to) = 6{x° - x + (t - t )u(x,t)). (35) 

It satisfies the adequate inviscid transport equation: 
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+ ^(u(x + ,t) + u(x ,t)j d x 



V re tXx,t\x°,to) 



(36) 



with the appropriate boundary condition: 

"P r et.OMkVo)|t=to =S(X-X ). 

Eq.(p6]) is valid inside correlation functions. Note that V re t. does not satisfy the naive transport equation 
(H) with k = 0, since eq.(^l|) yields: 

d t + u(x,t)d x ] Vret.(x,t\x°,t ) Si -{t-tof e(x)5"{x a -x+ (t-ta)u(x,t)) £ 

where the left hand side does not vanish due to the dissipative anomaly. To prove eq. (|36|) , let us expand 
V re t.( x iM x ° >*o) m Fourier series as: 



V re t.(x,t\x ,t ) = 
Plugging Vk into eq.(|36"|) gives: 



dk 
2^ 



r ° V k (x,t) with V k (x,t) = e -*(*-(*-*oM*.*)) 



d t + 



x+)+u(x-))d x 



V k (x,t) = (ik) 

+ (ik(t - 1 )) 



u(x) — —{u{x + ) + u(x )) 



Vk(x,t) 



d t u{x) + ~{u(x + ) + u(x ))(d x u(x)) 



Vk(x,t) 



The first term in the r.h.s. vanishes since correlation functions of the velocity field without derivative are 
continuous whereas the second vanishes thanks to the equation of motion (Eq). 
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• The forward probability distribution is: 

Vadv.(x,t\x°,t ) = -d x H(x,t\x°,t ) 

= (l - (t-t )d x u(x,tfj S(x° -x + (t-t )u(x,t)) (37) 

It satisfies the transport equation: 

V ad v.{x,t\x° ,t ) £* (38) 

which corresponds to the limit k — + of eq.(|6|). Remark that the Jacobian (1 — (t — to)d x ii(x,t)) is always 
positive since away from shocks d x u = X/t < 1/to and that on shocks d x u is negative. 

Eq. ( |3~7j ) implies that the forward probability distribution may be decomposed as the sum of the backward 
probability distribution plus an operator which is localized on shocks. Namely: 

V adv .(x,t\x°,t ) = jrret.(x,t\x ,to)-V(x,t\x ,t ) (39) 

with T>(x,t\x° ,t Q ) = (t — to) (d x u — |) V re t.(x, t\x°, t ) whose shock representation is: 

V(x,t\x°,t ) = ^2 * (f^ - l& ~ ^^ol) 5 ( x ~ x s(t)) with v x , X0 = - ^° 

with x(^) the characteristic function of the constraint C. Here the constraint may also be written as 
U J < v x,x < u t which means that the speed of the trajectory going straight from (xo,to) to (x,t) is 
between the two extreme values of the velocity at the shock. 



+ d x ^u(x + ,t) + u(x ,t)j 



6 Lagrangian trajectory distribution functions. 

The purpose of this section is to derive explicit formulae for the advanced and retarded one and two point 
function distributions of Lagrangian trajectories. We use these results to compute the short distance be- 
haviour of these correlation functions, the probability that a particle meets a shock or that two particles get 
glued together. We conclude with remarks on persistence problems. 



6.1 One-point functions. 

The one-point probability law of the velocity field u 

exp 



u(x, t) is: 



du 



This is well known since Kida but rederived for completeness in appendix A. Thus the one-point p.d.f. 
for backward and forward trajectories coincide and are equal to: 



Pl 1 l(x,t\x a ,to)=Pl 1 lix,t\x a ,t ) 



[i] 



t 



exp 



t(x-x°f 



2(<-*o) 2 



(40) 



It simply reflects the diffusion of the trajectories with ((a; — a; ) 2 ) ~ (t — to) 2 /t. For large t/to, this is just 
the ordinary dispersion of Brownian motion. But when t — to is small compared to to, the dispersion grows 
linearly with time because with high probability no shock has been met. 

It is instructive to compare this to the probability distribution for a particle starting at xo at time to to 
flow to x at time t without hitting any shock. As computed in Appendix B, this is equal to: 



t 



Mt-to) 2 



exp 



t(x - Xp) 2 
2(t-to) 2 



dx 



(41) 



In particular the probability that a particle does not meet a shock between to and t is to It. 

The probability distribution for a particle starting at xq at time to to flow to x at time t hitting exactly 
n shocks is more complicated for n > 0, and it is funny that the resummations for all values of n leads to 
such a simple result. 
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6.2 Two-point functions. 



The two-point trajectory p.d.f 's are slightly more lengthy to compute. The two-point p.d.f. for the velocity 
field Ui = u(xi,t) and u 2 = u[x 2 ,t) are recalled in appendix [a|. In the sequel, F t (z) stands for a variant of 
the error function defined by: 

F t (z) = e 2t / e 2t du 

J — oo 



• Let us first look at the backward probability distribution. Recall that it may be computed by a simple 
change of variables from the velocity distribution function. Thus for x\ > x 2 : 



P™ix,t\x ,t ) 



i 2 



{t - t ) 2 

At 
(t - to) 2 



S(A -t(vi- v 2 )) 



6(A - t{ Vl - v 2 )) 



1 



F t {-tv 2 ) + F t (tV!) 



dz 



(42) 



Note that as expected PjJ vanishes for t{x1 — x 2 ) > ^(^l ~ x 2) for (x\ — x 2 ) > 0. See comments below. 
Note also that in the coinciding limit X\ = x 2 one has: 



pW t Xx,t\x°,t )\ Xl -_ 



S{X ° - 4) y 2n(t-toy 6XP 
6(xl-x° 2 ) P^ t Xx,t\x°,t ) 



t(x-x ) 2 
2(t-i ) 2 



This means that two trajectories at identical final positions did start at identical initial points. The same 
applies to the n-trajectories probability distribution functions: 



Pjl (^14*0)1*,,=^ = 8(xl-xl^)Pl n e ; 1] (x 3 A^ v to) 
In other words, Lagrangian trajectories are statistically well-defined backwards. 



(43) 



[21 

• Consider now the forward probability distribution PIL. ■ It is less straightforward to compute, but 
the relevant information can be extracted from the formula for H^(xj, t\x®, to), which is a sum of two 
contributions: 

H PI ( Xj , 1 1 x° 3 , t ) = K ! ( X j , 1 1 x° 3 , t ) + K 2 [ Xj , 1 1 , t o ) 
Since iJPl is symmetric, it is enough to evaluate it for X\ > x 2 . To simplify the notations, we set 



Xj — x" 

A = Xi — x 2 > and = v„ . _o = 

3 x J ,x j t _ to 



Then 



and 



Ki(xj,t\x°,t ) = 



dz 



<tVl - A,t„ 2 +4) F t (f + z) + F t (f - z) 



if 2 (^,*N?)*o) = -r / / dz 2 6(A - (zi - z 2 )) 



I.V2 



Z2 + -. 



dz- 



1/2 2 \ k A 



a [F t (f + z) + F t (f -z)] 2 



(44) 



(45) 



Of course, one could recover the results for the backward probabilities using eq.([33|). The explicit use of 
eq.(p4[) leads to formulae for the forward probability which are not really illuminating. However, if PI can 
be interpreted as the probability that two particles starting at time to at points x\ and x\ respectively have 
abscissae at t larger than x\ and x 2 respectively. And indeed, one can check explicitly on the above formula 
for Jf PI many expected physical properties of trajectories: 
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- Particles do not cross each other: if < & (and in particular if (x® — x 2 )(xi — x 2 ) < 0), 

reduces to a one particle distribution: 



HM{x h t\x%t ) = 



' I due- 2 *' 2 

max(wi ,112) 

H^{xi,t\x%t ) for xi > x 2 
HW(x2,t\x 2 ,t ) for Xl <x 2 



(46) 



Taking derivatives with respect to x\ and x 2 , one finds a vanishing probability density if the respective orders 
of the particle positions have changed between initial and final times. 

- Trajectories are well defined forward: for fixed x\,x 2 and t, formula ( fl6| ) is valid for \x® — x 2 \ small 
enough, and leads to 



lim H [2 \x u x 2 ,t\xl,xl,t ) 



if[ 1 l(max(xi, x 2 ), t\x°, to) 

due~ uH ' 2 . 



*(si,s 2 )- 



Taking the derivatives with respect to X\ and x 2 gives: 

,[2] 



lim o P l2 i(x 1 ,x 2 ,t\x° 1 ,xlt ) = S^-x^lXx^M) 



8(xi - x 2 ) 



27r(i - t f 



exp 



t(x- x°) 
~2(t-t y 



0\2T 



(47) 



• Contrary to the backward p.d.f., P}$ v is singular at coinciding points: assuming x\ ^ x 2 a direct 
computation shows that 



P^ v Xxj,t\x° j ,t )=R{x 1 ,t\aP j ,to) 6{ Xl -x 2 ) + 



(48) 



The dots refer to terms regular at X\ — x 2 . The coefficient R(x, t\x®, to), which has dimension of the inverse 
of a length, is the probability density of aggregation of trajectories at point x. It is equal to: 



R(x,t\x1,x° 2 ,t Q ) 



2nt 



Ftitvx) + F t (-tv 2 ) 



for v\ < v 2 



(49) 



Let us note that this gives also the probability that n particles have collapsed, if v\ and v 2 refer to the speeds 
of the particles with the extreme initial positions. 

This formula simplifies if one is simply interested in the probability that two particles starting at distinct 
points at to have glued together at time t : integration over the final position gives for the total gluing 
probability 



t-t 



dv 



V/4 



(50) 



For fixed to and x\ — x 2 , and t — > 00 , this behaves like 1 



I*? 



x 2 \/y/Trt — to/t. This shows that distinct 



particles are sure to be at the same point at a late enough moment of the evolution, but this gluing occurs 
rather slowly. 

Another special case where the general formula simplifies is the limit of identical initial positions x\ — 



..0 _ »o 



x leading to 



R\x«= x o = V exp 



0\2 1 



t{x - x") 
'2{t-toY 



This should be compared with formula (47). The difference is exactly equal to the probability to go 
without shock from x° to x in the time interval [to,t] (fid]), and this has a good explanation: two particles 
starting at the same point stay sticked together, two particles starting at distinct point may coalesce only 
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when they meet a shock, so the difference in coliapse between starting at the same point and starting infinitely 
close is simply encoded in the probability that a single particle has met no shock. Integration over the final 
points shows that the probability for two infinitely close particle at time to to have glued together at time t 
is 1 — to/t. 

The collapse probability R may be computed in other ways. One way consists in using the operator 
product algebra of operator localized on shocks, cf. eq. (|2^) . Indeed, in view of the decomposition ( j39|) of 
T- > adv.( x >t\ x ° \to), one has the following operator product expansion: 

r Padv.(x 1 \x°)Vadv.{x 2 \x2) = T>(x x \x°) V(x 2 \x%) + regular 

= nfalx?) 5(xi - x 2 ) + ■ ■ ■ 

with 1Z the operator localized on shocks whose form factor is the product of those of the operator T>(x\x°), 
i.e.: 



K{x\x°) = X ( % > max |£, - vA S(x - x s (t)) 



2t 

Clearly, usinf gthe shock distribution recalled in the appendix, one gets (TZ(x\x^)) = R(x , t\x® , to) as com- 



puted in eq.(49) 



Another way to compute R(x,t\xj,to) is as follows. We know that particles do not cross, so that from 
the equation of trajectories, we can infer that two particles starting at x® and x 2 (x® > x%) respectively are 
glued together between x and x + dx at time t if and only if x — x 2 — (t — to)u(x,t) < (i.e. u(x, t) > v 2 ) 
and x + dx — x\ — (t — to)u{x + dx, t) > (i.e. u(x + dx, t) < vi). But as recalled in the appendix on shock 
distribution functions, the probability that u{x, t) > v 2 and u(x + dx, t) < v\ for v 2 > v\ is 



dx 
2^ 



dv+ dv.t(v + -v-)9{v+-V-)e- t{v + +v - )/2 - (51) 

J — oo 

This leads again to the above formula for R(x,t\x°,to)- 



6.3 A comment on persistence problems. 

To every random velocity distribution, one can associate domains on the x axis, defined as the intervals 
where the velocity v is continuous. Those domains change as shocks move and annihilate into other shocks. 
This is a typical situation where persistence concepts are useful. We have computed above two quantities 
that relate naturally to persistence. For instance the probability to move on a Lagrangian trajectory in the 
time interval [to, t] without meeting a shock, i.e. remaining in the same domain was found to be to/t. In the 
same vein, the probability for two particles starting on Lagrangian trajectories at distance x > from each 
other at time to to be at distinct positions at time t was found to be 

1-*^*° /" ^e- 2 /4, (52) 
t I icVt Wtt 

t-t 



which behaves for large t and fixed x and to as 

x to 
VttI t ' 



(53) 



In particular there is no unexpected persistence exponent. 

Let us note that the more usual definition of persistence, which is not related to Lagrangian trajectories 
but deals with points that do not move with time, leads to a different kind of behavior^] that can be computed 
by direct use of the distribution of velocities. For instance, the probability that a fixed point (say, the origin) 
is not hit by any shock in the interval [to,t] is 



/dy I dy' exp sup 
J \t'e[tn, 



2 /2 

y - y 



t'E[t ,t] 2t ' 



(54) 



°This comparison was suggested to us by Claude Godreche. 
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In the limit t/to — > oo, this exhibits the slightly nontrivial behavior 



• (55) 

quite different of the previous results for moving particles. 

7 Inverse cascade. 

We now consider properties of a tracer advected in the inviscid Burgers decaying turbulence. In particular, 
we argue that there is no dissipative anomaly and that the energy cascade is inverse. 

• As previously explained, in the inviscid limit the appropriate transport equations are eqs. (§,§). In the 
limit k — ► 0, their solutions may be written in terms of the backward and forward probability distributions. 
Namely, 

T(x,t) = J dx° V ret .{x,t\x°M) T (x°) (56) 
= T Q (x - (t - t )u(x,t)) 

and 

p(x,t) = J dx° V adv .(x,t\x ,t ) p (x°) (57) 
= (!-(*- t )d x u(x, t)) p (x - (t - t )u(x, t)) 

where To (a; ) and po(x°) are the initial conditions at time to. 

Since correlations of the trajectory probability distributions are computable, there are not much difficul- 
ties to evaluate correlations of the scalars. Let us illustrate this by showing that there is no dissipation of 
energy for the tracer T(x, t) and hence no dissipation anomaly for T. The mechanism for that property is 
similar to the one described in ref. in the case of the compressible Kraichnan's model. Assume that one 
is given the translation invariant two-point function of the initial data: 

{T (x 1 )T (x 2 )) = T(xi - x 2 ) 
The density of energy of the tracer is £(x, t) — ^T 2 (x, t). Its average is: 

(£(t)) = lJdx°dxlP^i(x,x\xlxl)(T (x°)T (x^ 
= \J dx a Pl 1 i(x\x°)T(0) = ^T(0) 

where we have used eq. (fl3|) for P^ t at coinciding points and the normalization condition for . Thus 
energy is conserved in mean, (£(t)) = £q, and this is due to the fact that the trajectories are statistically 
well-defined backward. Notice however that at fixed initial data the density of energy decreases at large time 
as (£(x,t)) ~ J dyT^y) if the integral converges. 

More generally, the well-defined character of the trajectories may also be formulated as the following 
operator product identity: 

Vret.{x,t\x\,t Q )V ret .{x,t\Xx,-to) = tf(aff? - 3?§) Vret. {x, t\x\, t ) 

As a consequence, any products of solutions of the transport equation (^) at k — will also be solution. In 
particular, any powers of T(x, t) are also solutions: 



d t T n (x, t) + ~ (u{x+, t) + u{x-,tj)d x T n (x, t)=Q 



inside correlation functions. This shows the absence of dissipative anomalies in the passive free advection 
which means that the fields nT n &iT vanish inside correlation functions at k — 0. 
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• This is the sign of the absence of a direct energy cascade, as in the two dimensional turbulence in which 
the energy cascade is inverse, i.e. toward the large scales To show it more explicitly let us now assume 
that one is injecting energy to the tracer such that the transport equation is now: 

d t T{x,t) + ~(u(x + ,t)+u(x-,t))d x T(x,t)=f(x,t) (58) 
with fix, t) the forcing term. Solutions of this equation with zero initial data at time to are: 

T(x,t) = f ds [ dyV ret ,(x,t\y,s) f(y,s) (59) 



t J 
t 

ds fix — (t — s)u(x. t), s) 
Assume that the two-point function of the force is delta-correlated in time: 

(f(Vl, Sl)f(V2, S2)) = C L (yi ~ V2) S(S1 - 82 ) (60) 

with Cl{x) a smooth function varying on scale L and with rapid decrease at infinity. The energy injection 
rate is e = |Cx(0). Using again the fact that trajectories are well-defined backward, eq. (ff3|) , one finds that 
the average of the tracer energy density at time t is: 



(S(t)) = IJ^dsJdy P^{x,t\y,s) C L (0) 



'to 

= \(t-ta)C L (0) = (t-to)e (61) 

where, again, we used the well-definedness of the trajectories (see eq. (fl3|)) and the normalization of the 
probabilities. Thus the total amount of energy injected into the system is transfered without dissipation. 

To decipher in which mode the energy is injected, let us consider the scalar two-point function at distinct 
points. For forcing delta-correlated in time as in eq.(pOf), the two-point function is: 



(T(x 2 ,t)T(x 1 ,t)) = / dsdy 1 dy 2 Pl 2 } t .{x j ,t\y j ,s) C L {yi- y 2 ) 
It behaves at large time and fixed positions as: 

rl ds 

(T(x2,t)T( Xl ,t)) ' ' ' ' 



f ds 

ty/n - \x 2 - Xi\Vi ) / — = C L (s\x 2 - xi\ 

' JO V7T 

+ F(x 2l xi) + 0(1/V*) 



with F{x 2 ,x\) finite as t — > oo scaling as \x 2 — x\\ at small distance. The energy is thus transfered to the 
mode corresponding to the first line of the above equation. Its amplitude increases with time. It is a soft, 
although non constant, mode varying smoothly and slowly. 

To make manifest the absence of dissipation, consider products of the forced scalar ([39]) at coincident 
points. One has: 



T n (x, t) = l[ / ds 3 f{x - (f - 8j )u(x, t), Sj ) 

3=1 Jt0 



Using again eq.(|36|) or (28), one deduces that inside correlation functions: 

d t T n (x, t) + i (u[x+ 1 1) + u(x~,t)) d x T n (x, t)=n fix, t) T"" 1 ^, t) 

This shows that there is no dissipativc anomalies at k = in the scalar advection. Note that what we have 
described is a limit when k goes to zero first and then t goes to infinity to reach the stationnary state. 

Again the mechanism is similar to the one found in compressible Kraichnan's models the injected 
energy is accumulated in the soft mode, there is no dissipative anomaly and the energy cascade is inverse. 
This is directly related to the fact that the trajectories are statistically well-defined. 
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A Velocity probability distributions. 



In this appendix, we recall known formulae for the one point and two point probability distributions for 
velocities (see eg.|Bj). We just give a reminder of the computational rules and illustrate it in the case of the 
one point velocity p.d.f. A further illustration is given in appendix |B[ 

We define S(x,t) — miiij + ^ x ^ ^ so that u(x,t) ~ d x S(x,t). The pairs {4>j,yj) are described by 

a Poisson point process, saying that the cell of size dipdy in the (0,y)-plane is occupied with probability 
'dtfidy , disjoint cells being independent. This leads to the following useful fact that if D is any measurable 

s e. i„ tire ,,-piane, the probability that aii celis in D are e m p. y is /. We cal, that pair («,,„) 

giving the minimum of S at the point (x,t) the parameters at (x,t). 

A.l One point velocity p.d.f. 

We look for the probability P(u(x, t) S [v, v + dv]). The law for the Poisson point process reads in this case : 
— A cell ((/>, y) with € [v 7 v + dv] is occupied. 

- The cells in D = {(</>', y') such that 0' + ( *~*f )3 < cj) + = + 2^} are empty, 

so 

P(u(x,t) e[v,v + dv]) = [ e*d<j)dye~ $o e * d4, ' dy ' . 

Let us do this computation in detail. First we do the integral over 0', which varies between — oo and 
+ i£ ir l! - This yields 

C j C „&+(x-y) 2 /2t-(x-y' ) 2 /2t 

P(u(x,t) € [v,v + dv]) = eUcf>dye-J e dy . 

J ^e[«,u+(i«] 

Then, we integrate over </> to get 

/" p -( a; -y) 2 /2t 

P(u(x,t) € [v,v + dv]) = / 



^g^+cfe] y Je-(-y') 2 /2t dy r 

Let us note that the possibility to integrate explicitly over the variable 4>' parametrizing the empty domain 
D and over the "center of mass" of the variables 4> parametrizing the occupied cells is typical. In this explicit 
example, the other integrations are also immediate, but this is rather unusual. 

The y' integral gives a factor l/-\/2nt and the integration domain for y is infinitesimal, so y — x — vt and 
dy = tdv. Finally: 

P{u{x,t) e [v,v + dv]) = \U-e~ tv2 / 2 dv. 

V 2tt 

Let us observe that it has total mass 1, ensuring that this computation, which does not take shocks into 
account, does not miss any event of nonzero measure. This is a sign that shocks are diluted. 

A. 2 Two point velocity p.d.f. 

We look for the probability P(u(xi, t) G [v±, v\ + dvi], u(x2,t) € [i>2, «2 + dv2\)- By symmetry, we may (and 
shall) assume x\ — X2 = A > 0. There are two possibilities : 
(*) one parabola: 

— A cell (4>, y) with Xl ~ v g [vi,vi + dvi] and X2 ~ v G [v2,V2 + dv2] is occupied. 

- The cells in D = {(</>', y') such that </>' + (xi ^' )2 < 4>+ or </>' + (x ^ f < 0+ to^} are empty. 
(**) two parabolae: 

— A cell (4>i,yi) with Xl ~ yi <E [v\, V\ + dv\] and a cell (</>2, Z/2) with X2 ~ V2 g [v 2 , v 2 + dv 2 ] are occupied, such 
that fa + <fa + and fa + < fa + 

- The cells in D = {(0', y') such that fa + ^'f )2 < fa + {xi ^ l)2 or fa + ( " 2 ^' )2 < fa + ( ' 2 ~f )2 } are 
empty. 

Accordingly, P(u(x\,t) e [v\,v\ + dv\], u(x 2 , t) S [v2,v 2 + dv2]) is a sum of two contributions Pi pa rab and 
-P2parab which are found after some computation to be: 

Piparab = t 2 dv 1 dv 2 8{& - t(v x - V 2 ))—, , 1 , , 

F t {-v 2 t) + F t {vit) 
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and 



rtv 2 +A/2 z 2 /t 
P 2 parab = Atd Vl dv 2 6(A - t( Vl - W2 )) e -*(^I)/2+A 2 /4t / dz 

Jtvi-A/2 



F t ($-z)+F t ($+z) 



Let us recall that F(x) = f*^ dye~ y '' ''/ 2 and F t (x) = s/2rtF (x / . 

Again, one can check explicitly that the sum has total mass 1, or even better that the integral over v\ or 
V2 gives again the one point p.d.f. This computation shows that Pi pa rab which lives on a codimension one 
hyperplane, is completely determined as a kind of boundary of P2parab- 

A. 3 Distribution of shocks. 

The two-point p.d.f for velocities allows to compute the probability to have a shock such that u(x, t) = v + 
and u(x + dx, t) — v- between x and x + dx by taking A — > 0. The result is 

^dv + dv_ t{v+ - v_)6(v+ - v^e'^+v-)/ 2 - (62) 

ZTT 

This can be expressed as the probability to find a shock of amplitude fi/t = v+ — v- and velocity £ = 
(v+ +V-)/2 in the interval [x,x + dx] as 

^«^(^)e-« 2t ^ 2/4t - (63) 

In particular, the probability to have a shock in the interval \x,x + dx] is dx/sfwi. This involves only 
configurations with two parabolae whereas the probability that there is no shock in a finite interval [x, x'] is 
computed with configurations involving one parabola and found to be J dy(F t (y — x) + F t (x' — y)) . This 
makes it intuitively (if not mathematically) clear that with probability one a finite interval contains only a 
finite number of shocks. 



B One-point p.d.f. without shock. 

In this appendix, we compute the probability P(x, t\x , t ) no shock that a particle starting at point x at 
time to arrives in [x, x + dx] at time t without ever meeting a shock. This corresponds to the following 
configuration: 

- A cell {4>, y) with x + x ° t ~ v (t — t ) € [x,x + dx] is occupied. Let v = (x — y)/to- 

- The cells in D - {{<f>', y') such that <f>' + -J°)-v'f < + -t a )- y ? for somc ^ g &re 
empty. 

The second constraint seems complicated. We claim that it is equivalent to the extreme constraint for 

t' = t: 

- The cells in D = {(</>', y') such that 0' + (*°+^-M^') 2 < ^ + {x +v{t-t )-yf y &ic cmpty _ 

This is a direct consequence of an important property of trajectories. As already stated before, Lagrangian 
trajectories stick to shocks as soon as they meet one. We can even be a bit more precise. Suppose that at 
time (x,t) the parabola of parameters {4>,y) dominates (<f>',y'), so 

^+^^<0+^^, (64) 

or better 

(f> _ (f)l< (y-y')(2x-y- y >) (65) 

Consider a fictive particle moving at constant speed v = (x — y)/t and arriving at point x at time t. At time 
to < t it was at point Xo = x — v(t — t ). The identity 

{y - y'){2x - y - y') (y - y')(2x - y - y') _ t - tp ,, 2 , , 

2t 2t 2tt Q yu g > y > 
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proves that the parabola of parameters ((/>, y) was already dominant at (xo,to). This proves that the equiv- 
alence of the two above definitions of the forbidden domain D. This means also that if point x is not on 
a shock at time t and u{x,t) = v, there is a unique backward Lagrangian trajectory through (a;,t), defined 
back to time t$ and such that at to the particle was at point xq = x — v(t — to). 
So we need to compute 

Jx + ^2-(t-to)£[x,x+dx] 

Again, integration over <f>' ', <fi and y' is straightforward, and yields 

P : . (*, t\x , t Q ) dx = (I) ^/ 27r(t l to)2 
Hence, the probability that no shock is met in the interval [t ,t\ is simply to/t. 



t(x - x y 
'2(t-to) a 



dx 
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